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■ 1, Abstract. Optical images and high-dispersion spectra have been obtained of the ejected material surrounding the 

' pulsating AGB star Mira A. The two streams of knots on either side of the star, found in far ultra-viollet (FUV) 

' GALEX images, have now been imaged clearly in the light of Ha. Spatially resolved profiles of the same line 

• , reveal that the bulk of these knots form a bi-polar outflow with radial velocity extremes of ± 150 km s~^ with 

■ respect to the central star. The South stream is approaching and the North stream receding from the observer. A 

displacement away from Mira A between the position of one of the South stream knots in the new Ha image and 
^ , its position in the previous Palomar Observatory Sky Survey (POSS I) red plate has been noted. If interpreted 

^ ' as a consequence of expansion proper motions the bipolar outflow is tilted at 69° ± 2° to the plane of the sky, 

C/) ' has an outflow velocity of 160 it 10 km s~^ and is ^ 1000 y old. 
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1. Introduction 

The binary system Mira AB has attracted attention over 
400 yr (for a summary see Hoffleit, 1997), for apart from 
being very close to the Sun (the Hipparcos distance is 
107 ± 10 pc) Mira A is the archetypical pulsating vari- 
able star with a period of hundreds of days and brightness 
range of 8 mag (Reid & Goldston 2002). Its periodic out- 
bursts over the last 80 cycles are listed by Percy & Bagby 
(1999). Hubble Space telescope images have shown Mira 
B to be separated by ^ O.b" from A (= 54 AU; Matthews 
& Karovska 2006). Mira A of intermediate mass is cer- 
tainly in its volatile asymptotic giant branch (AGB) phase 
whereas B was previously identified as a White Dwarf but 
recently as a somewhat innocuous F-type star similar to 
the Sun though with an accretion disk formed by the mass 
ejections from A (Karkova et al. 2005; Ireland et al. 2007 
and refs. cited therein). In any case B is too well separated 
from A to affect significantly the overall dynamics of the 
ejected material. 

Interest in Mira A has been further heightened by the 
discovery with the GALEX satellite of a 2° long 4 
pc) FUV emitting tail (Martin et al. 2007). They propose 
that this traces the interactions with the local interstel- 
lar medium (ISM) over the past 30,000 y of the ejected 
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material from Mira A. This is plausible (Wareing et al. 
2007) for the star has a high space velocity of 120 km 
s~^ at an angle of 28° with respect to the plane of the sky 
tilted away from the observer and what appears to be a 
preceding bow-shock (Raga & Canto 2008). Perhaps rele- 
vant is the tail of similar dimensions apparently projecting 
from the similarly mass-ejecting proto-typical Luminous 
Blue Variable (LBV) star P Cygni (Meaburn et al. 1999); 
though later, deeper, observations in the higher excitation 
[O III] 5007 A line show this tail as part of a larger complex 
of filaments but still convincingly associated with the star 
(Meaburn et al. 2004; Boumis et al. 2006 and references 
therein). 



In the GALEX FUV images a group of emission knots 
close to Mira AB and called North Stream and South 
Stream by Martin et al. (2007) have been identified as re- 
cent manifestations of the mass ejections from Mira A. As 
these are within a field of a few arcmin diameter they have 
now been investigated at optical wavelengths using occult- 
ing strip imagery, to exclude scattered continuum light 
from the bright stellar image, and longslit, spatially re- 
solved, high-dispersion spectrometry. The results of these 
optical observations will now be presented. 



2 



J. Meaburn et al.: High speed outflows from Mira 



2. Observations and results 

Observations were made with the Manchester Echelle 
Spectrometer (MES-SPM - see Meaburn et al. 1984; 2003) 
combined with the 2.1-m San Pedro Martir, (B.C. Mexico) 
telescope on 2008, 26 to 28 Oct. A SITe CCD was the de- 
tector with 1024x1024, 24 fim pixels although 2x2 bin- 
ning was employed. The 'seeing' throughout was ^ V . 

2.1. Imagery 

MES-SPM has a limited imaging capability with a re- 
tractable plane mirror isolating the echelle grating and a 
clear area (4.37x5.32 arcmin^) replacing the spectrometer 
slit. 

The image in Fig. 1 is a subset from this larger field ob- 
tained on the 2008 28 Oct. An occulting strip of chromium 
suppressed the bright image of Mira. The integration time 
was 2000 s. The coordinates (J2000) were added using the 
STARLINK GAIA software. The 90A bandwidth interfer- 
ence filter transmits predominantly the Ho^ nebular line as 
[N II] 6548 & 6584 A emission has been shown (Martin et 
al. 2007) to be very low. Confusing, faint star images were 
ehminated using the PATCH routine of the STARLINK 
GAIA software. Deep and lighter representations of the new 
Ha image are compared in Fig. 1 to the FUV GALEX im- 
age from Martin et al. (2007). 

2.2. Longslit spectroscopy 

Spatially resolved, longslit Ha line profiles were obtained 
with the MES-SPM along the lines marked 1-6 in Fig. 
2. These are only partial lengths of the full NS slits as 
measurable Ha emission only occurred over small sections 
of the full slit lengths containing the bright knots in the 
North and South streams. The increments along the slit 
length each corresponds to 0.63''. 

In this spectroscopic mode MES-SPM has no cross- 
dispersion consequently, for the present observations, a 
filter of 90 A bandwidth was used to isolate the 87*^ echelle 
order containing the Ha and [N ll] 6548 & 6584 A nebular 
emission lines. The slit width was always 300 /im which 
is = 3.8'' on the sky and 20 km s~-^ spectral halfwidth 
(HPBW). Each integration time was 1800 s. Also indicated 
in Fig. 2 are the directions of the 2° long FUV emitting 
tail detected by Martin et al. (2007) and of the proposed 
bowshock found in the same wavelength domain. 

The longslit spectra were cleaned of cosmic rays, bias 
corrected and calibrated in wavelength to dz 1 km s"""^ 
accuracy in the usual way using STARLINK FIGARO soft- 
ware. The greyscale representation of the position-velocity 
(pv) arrays of Ha line profiles (after subtraction of the sky 
background spectra) for the partial slit lengths shown in 
Fig. 2 for Slits 1-4 and 5-6 are shown in Figs. 3 & 4 re- 
spectively. Here contours of the Ha relative surface bright- 
nesses with linear intervals have been overlain on the neg- 
ative greyscale representations. As no standard star was 
observed the absolute surface brightnesses are unreliable 



and will not be presented here. Line profiles from vari- 
ous incremental lengths of these pv arrays of profiles are 
shown in Fig. 5. 



3. Discussion 

3.1. The bi-polar outflow 

The radial velocities of the South stream of knots shown 
in Fig. 3 show a shift of around -150 km s~^ with respect 
to the heliocentric systemic radial velocity V^ys = 56 km 
s"^ as given by Matthews et al. (2008 - N.B. Vlsr = 
Vrel - 9-44 km s~^ for Mira A). A very similar shift, but 
now to positive radial velocities can be seen in Fig. 4 for 
the North Stream on the opposite side of Mira A. Only 
the more distant arc from Mira A (Figs. 1 & 2) covered by 
the element Slit 5a (see the profile in Fig. 5) has a greater 
shift (200 km s~-^) in radial velocity with respect to Vsys- 
There seems therefore to be a strong bi-polar structure to 
the principal elements of this outfiow. 

This expansion radial velocity measurement of the bi- 
polar lobes has to be combined with a measurement of 
their expansion tangential velocity relative to Mira A to 
determine their actual expansion velocity. The expansion 
proper motions (PMs) then of the knots relative to Mira A 
are required for this purpose. Simple inspection of the Ho; 
image in Fig. lb and FUV image in Fig. Ic reveal that the 
knots in Ha emission are ~ 2" closer to the reference star 
at the bottom of both frames than in the FUV emission. 
This could indicate the displacement caused by a high 
expansion PM in the intervening 1.9 y between the two 
observation leading to an expansion velocity of ^ 600 km 
s~-^. However, a more likely cause is the spatial diff'erence 
within the knots of the Ha and FUV emission regions. 

This suggestion is born out by the spatial displacement 
of Knot B (Fig. 2) in the first Palomar Observatory Sky 
Survey POSS I red plate taken on the 26 Nov. 1954 and 
the Ha image in Fig. lb obtained on the 20 Oct 2008. 
Only Knot B is clearly identifiable on this POSSI plate 
with the image of the brighter Knot A confused with that 
of a faint star. Curiously none of the Mira knots nor this 
star are detectable on the follow-up POSS II (23 Aug 
1995) red plate. However, using the STARLINK GAIA soft- 
ware the displacement of Knot B, measured with respect 
to the faint star field, in the intervening 53.90 y is +2" 
in Right Ascension (a) and -17"in Declination {S) with 
an uncertainty of dz2" in both measurements determined 
mainly by the poor quality of the image of Knot B in the 
POSS I plate. This low quality is not surprising for the 
POSS bandwidth is many times greater than the present 
one, photographic emulsion was the detector and the fo- 
cal length shorter. The PM of Knot B is then measured 
as d(a,^)/dt = (+37, -319) mas y"-*^ each to around ±40 
mas y"-"^ uncertainty. Here it is assumed that Ha domi- 
nates in the POSS I detection even though the bandwidth 
also encompasses the [S ll] 6717 & 6731 A lines and much 
more continuum light than the present observations. 
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For comparison a similar PM measurement of Mira A 
itself was made between the POSS I and POSS II red 
plates. In the intervening 40.74 y the displacement is - 
1.1" in right ascension and -8.4'^ in declination to give 
d{a,5)/di = (-40, -206) mas y"^ with (±30) mas y-^ un- 
certainties for the Mira A positions for these were only 
determined from the centroids of the diffraction spikes in 
the two images but still corrected for the small diflFerences 
in coordinates of the faint star field. 

It is firstly interesting that from these PM measure- 
ments Knot B is moving along PA = 173° which is reason- 
ably aligned with bipolar axis of the South stream knots 
in Fig. 1 and tilted by a few degrees to the direction of 
Mira A (as noted by Martin et al. 2007) along PA ^ 190° 
which is also the orientation of the long FUV tail (Fig. 
2). An expansion along this bi-polar axis of ^ 115 mas 
y~^ of Knot B relative to Mira A is therefore indicated to 
give an expansion tangential velocity of 58 ± 20 km s~^. 
When combined with the expansion radial velocity of 150 
km s~-^, this indicates an outflow velocity of ?^160 km 
s~^ tilted for the southern lobe at 69 ± 15° towards the 
observer with respect to the plane of the sky when all un- 
certainties are considered. More precise PM measurements 
of the knots are required to refine these estimations. The 
PM of Knot B with respect to Mira A, which is 127^' away 
from it, gives a time of ^ 1000 y since its ejection. The 
Mira A outburst in 1070 (Ho 1962; Hoffieit 1997) then 
become candidate for the creation of these bi-polar lobes. 

3.2. The leading bow-shock 

The preceding arc of FUV emission in the direction of 
motion of Mira A has been interpretated as a bow-shock 
being generated by the 120 km s~^ space motion of Mira 
A through the local ISM (Martin et al. 2007). It should 
now be considered that the bi-polar outflow identifled here 
could contribute to the creation of this bow-shock. In any 
case the radial velocity measurements (Figs 3 & 4) and 
the angular extent of the bi-polar lobes (Fig. 1) certainly 
indicate, without consideration of expansion PMs, that 
the outflow velocity is ~ 200 km s~^. Incidently sev- 
eral planetary nebulae (PNe) with central binary systems, 
one star of which has evolved beyond its Mira A AGB 
phase, exhibit bow-shocks emitting optical lines, as they 
plough through the local ISM (Bond & Livio 1990 and refs. 
therein). However, the question is posed that if bi-polar 
ejecta are favoured as the creators of these bow-shocks 
then why has GALEX only detected one such system? 

3.3. Evolution 

Mira A is currently in its AGB phase and shows all the 
signs that it could evolve soon into a complex PN maybe 
with multiple bi-polar lobes. For comparison, such lobes of 
the PN NGC 6302 have been conclusively shown to have 
been ejected at a high velocity (630 km s~^; Meaburn et 
al 2005 & 2008) and formed over a timescale of ^ 2000 



yr. The bi-polar outflow from Mira A revealed here could 
easily develop into such features over a similar timescale. 
Meanwhile, if Mira A itself evolves into a hot PN progen- 
itor as it sheds its outer layers in its current AGB phase 
during this period these high speed lobes will become ra- 
diatively ionized and emit a wide range of emission lines. 
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Fig. 1. Deep and lighter representations in panels a) & b) respectively are of the new Ha image of the knots obtained 
on 2008 28 Oct. and compared in panel c) with the FUV image of Martin et al. (2007) obtained on 2006 18 Nov. - 15 
Dec. The occulting strip seen in a) & b) minimizes the scattered light of Mira A within the telescope and instrumental 
optics. Coordinates throughout are 2000 epoch. 
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BOWSHOCK 

Fig. 2. Those parts of the lengths of the spectrometer slits containing useful Ha profiles are marked as 1-4 over the 
southern knots (South Stream) and 5-6 over the northern knots (North Sream). The directions of the FUV tail and 
bow-shock are also indicated. The top panel is of the new Ra image in Fig. 1 but printed lightly to show the details 
of the North Stream and the bottom window is deeper for the South Stream. The relative positions of the two panels 
in this figure are somewhat arbitary. 
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Fig. 3. Contoured and greyscale representations of the pv arrays of Ha line profiles from the South stream of knots 
shown in Figs.l & 2. These are for the slit lengths 1-4 shown in Fig. 2. The spectra of the background have been 
subtracted in all cases. The heliocentric systemic radial velocity Vsys of Mira A is arrowed. North is to the top of each 
slit length and each increment is = 0.63 " . 
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Fig. 4. As for Fig. 4 but for the northern knots covered by sht lengths 5 and 6. 
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Fig. 5. Sample line profiles extracted from the pv arrays of line profiles in Figs. 3 and 4 are shown. That for Slit la is 
for Ha profiles from incremental lengths 30 to 50 coadded along the Slit 1 in Fig. 3 and similarly the profiles Slit 4a 
and b are for lengths 68 to 83 and 31 to 37 respectively for Slit 4. The profile marked Slit 5a is for increments 132 to 
144 along Slit 5 in Fig. 4. As the background spectra have been subtracted in all of these profiles the noise fiuctuates 
around the zero level. Again the heliocentric Vsys is arrowed. The profiles have been smoothed by a Gaussian of 1.5 
times the 20 km s"""^ instrumental HPBW. 



